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Abstract. Recent progress in the understanding of the effect of electrostatics in 
soft matter is presented. A vast amount of materials contains ions ranging from 
the molecular scale (e.g., electrolyte) to the meso/macroscopic one (e.g., charged 
colloidal particles or polyelectrolytes). Their (micro)structure and physicochemical 
properties are especially dictated by the famous and redoubtable long-ranged 
Coulomb interaction. In particular theoretical and simulational aspects, including 
the experimental motivations, will be discussed. 
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1. Introduction 

Probably one of the most well known and understood ionic materials is sodium chloride 
(NaCl). In its solid form (i.e., NaCl cubic-like crystalline lattice), the experimentally 
measured heat of vaporization (7.92 eV) can be deduced (within about 10%) from a 
straightforward lattice sum of the form0 



leading to the theoretical Madelung energy (here E]\j = — 8.94eV) [Il[2]. This striking 
good agreement demonstrates that electrostatics is indeed the relevant ingredient 
governing our ionic crystal [3]. In its liquid form, NaCl plays a fundamental role in 
soft matter, since it controls the degree of screening of the Coulomb interaction in all 
water based solutions. It is exactly this type of problem that this review will address: 
Electrostatics in soft matter. 

Virtually all materials are more or less charged at the mesoscopic scale, depending 
on the degree of the polarizability of the embedding solvent (or matrix) and the solute 
particles (e. g., colloidal particles, polymers, membranes, etc.). The most well known 
example of polar solvent is evidently water which plays a crucial role in life, biological 
processes as well as industrial applications. When the solute particles are polar too, 
they can then dissociate into charged particles (also called macroions) and (microscopic) 
counterions. The counterion distribution near macroions turns out to be decisive for 
the surface properties of the latter. 

The pioneering works of Gouy and Chapman [H [5], realized almost one century 
ago, concern the counterion distribution near a planar charged interface. Applying 
the presently called Poisson-Boltzmann theory, they demonstrated that the counterion 
distribution profile decays algebraically as a function of the separation from the wall 
with a characteristic length that is inversely proportional to the surface charge density 
of the wall. Ten years later, Debye and Hiickel [6] accomplished a fundamental advance 
towards the understanding of screening. This theory originally developed for electrolytes 
(i.e. a solution of microscopic cations and anions such as Na"*" and Cl~) and based on 
the linearization of the Poisson-Boltzmann equation is now widely used in plasma and 
solid state physicsm 




Mean-field theories are appealing tools due to their intuitive and clear physical 
basis, and are robust theories as long as electrostatic correlations are not too important. 
In many practical situations (chromatin, polyelectrolyte multilayering, charged colloidal 

f The resulting energy in equation ([1]) corresponds to the cohesive energy per NaCl molecule. An ion 
(either Na+ or Cl~) is placed at the origin and aj = depending on the type of ion sitting at the 
lattice position rj. e = 1.602 x 10^ ^^C stands for the usual elementary charge, eo = 8.854 x 10~^^F/m 
for the vacuum permittivity, and a — 2.8lA for the NaCl lattice parameter. 

I Note that a similar potential of interaction (so-called Yukawa potential) arises at the subatomistic 
scale to describe the cohesion of the nuclear matter. Nonetheless, in nuclear physics, the interpretation 
of this potential in terms of screening is not adequate. 




(1) 
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suspension, etc.) electrostatic correlations are strong enough to make mean-field theories 
fail even on a qualitative level. Two striking and natural consequences of electrostatic 
correlations, that can not be explained by mean-field theories, are charge reversal (also 
called overcharging) and like charge attraction: (i) Overcharging concerns the situation 
where a macroion is locally covered by a cloud of counterions whose global charge 
overcompensates that of the macroion so that the net charge (or effective charge) changes 
sign; (ii) Like charge attraction is the counterintuitive effective attraction between two 
macroions carrying the same electric charge sign. 

A colloidal suspension, the classical material of soft matter science, can crystallize 
via a strong enough mutual electrostatic repulsion. An understanding of the 
resulting phase behavior necessitates approaches where particle-particle correlations 
must obviously be taken into account. This constitutes another example where 
approaches going beyond the mean-field level are required. 

The present work examines the problem of electrostatics in soft matter systems 
using simple theoretical models and computer simulations. The role of the little 
counterions is addressed in chapter [2l The relevance of excluded volume (i.e. the 
finite hard-core size of the constitutive ions) is discussed in chapter [3l The problem 
of image charges as occurring near curved dielectric interfaces is presented in chapter 
m The basic physics in more complex processes such as polyelectrolyte adsorption and 
multilayering is elucidated in chapter [5l Colloidal dispersions in strong confinement 
are presented in chapter [61 Finally, a conclusion and possible outlooks are provided in 
chapter [71 

2. Electrolyte at interfaces 

2.1. Foundations of electrostatic mean field theories in soft matter 

This part deals with the foundations of the electrostatic mean field theories in soft 
matter. It is written on a pedagogical level such that the non-specialist reader should 
be in a position to easily capture the underlying physics. Nonetheless, the expert will 
also certainly find some clarifying ideas in the forthcoming discussion. 

2.1.1. Poisson-Boltzmann theory The model system we have here in mind is sketched 
in figure [H We have to deal with a uniformly charged interface with a surface 
charge density a, separating the semi-infinite substrate from a simple electrolyte [i.e. 
univalent cations (-I-) and anions (— )] and whose bulk density is po- The system 
is globally electroneutral and the embedding solvent is merely characterized by its 
dielectric constant. In this context, the first theoretical determination of counterion 
distribution for an inhomogeneous fluid was realized by Gouy ^ and Chapman [5] 
independently almost one century ago. This mean-field approach corresponding to the 
so-called Poisson-Boltzmann theory is going to be now explained and discussed. 

A central quantity in the statistical mechanics of fluids is the potential of mean 
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z- distance from charged surface 



Figure 1. Model for a simple electrolyte near a (negatively) charged surface. 



force (PMF). The latter corresponds to the potential stemming from the effective force 
between two objects. The term "effective" means here a thermodynamical averaging 
whose form is dependent on the ensemble (e.g., canonical, grand canonical) under 
consideration. For the sake of simplicity we will consider the thermodynamical (i.e., 
macroscopic) limit where all ensembles are equivalent. 

A good starting point is provided by the exact Poisson equation relating the the 
mean electrostatic potential (MEP), tp{z), to the PMF Wa{z) as follows: 

Aij{z) = {exp [-(3w^{z)] - exp [-f3w^{z)]} , (2) 

where esoiv is the relative permittivity of the solvent (for water esoiv ~ 80), j3 = l/iksT) 
the reduced inverse temperature with ks being the Boltzmann constant and T the 
absolute temperature. The central approximation of the PB theory is to now set: 

w±{z) ±e'?/'(z), (3) 
such that the exact Poisson equation ([2]) becomes in the framework of the PB theory: 

Ai;{z) = ^ sinh [(3ei,{z)] (e = e^e^oi.) (4) 

which is the well known PB equation. The resulting MEP reads [7]: 

2kBT [1 + 7exp(-fi;2;) 

iIj[z) = In 

e [i- — 1 exp[—KZ) 

where 7 is given by the positive root of: 



(5) 



Y + (2k6)7 -1 = so that < 7 = -nb + V 1 + {nby < 1, (6) 

We have introduced here in equation and equation two important length scales, 
namely the screening length k~^: 

= SiriBPo (7) 
and the Gouy- Chapman length b 
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Figure 2. Reduced (Gouy-Chapman) counterion distribution /3+(z)(27r^B6^) = 
(i+z/b)" given by equation (fTT|) . It is precisely at z = 6, that the cumulated 
counterions (shadowed region) half-compensate the charge of the surface. In other 
words, the counterion integrated charge at 2; = is exactly —a/2. The strong Coulomb 
coupling limit (Moreira-Netz) p^{z){2ntBb'^) = exp(— z/6) as given by equation (fT9|) 
is also shown for direct comparison. 



where (.b is another third relevant length in charged soft matter known as the Bjerrum 
lengtfQ and reads: 

^« - sir- f'" 

The salt-free case can be actually easily obtained by considering k6 — > in equation 
([5]) and ([6]). Doing so we find: 

lim ^(2) = ^-ln[l + 2/6] + — ^In— . (10) 

The corresponding counterion distribution, p+(-2) = po exp [— /5e'?/'(2;)], is then merely 
given by: 

P+W = ^^J^. (salt -free) (11) 

It is this formula flTTl) that is referred to as the Gouy-Chapman counterion distribution. 
The corresponding plot can be found in figure [2l 

To better understand the physical meaning involved in the approximation ([3]), we 

f The physical interpretation of the Bjerrum length is straightforward: It is the distance between two 
elementary charges e that leads to an electrostatic interaction equating ksT . 
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shall make use of the exact so-called Yvon-Born-Green (YBG) hierarchy [S] that reads: 

ga(3iri2, Zi, Z2)pf3{z2)d^il2,) 



47reri2 



- ViWaiZi) = ^^^z - X] / ^1 

13=1 

Equation f[T21) can be seen as a "Newtonian" version of the statistical Poisson equation 
([2]). The left hand-side of equation (fT2!) represents the effective force felt by the test 
ion 1 of species a = ± at prescribed location ri = {xi,yi,Zi). The right hand side 
of equation f|T2|) is made up of two contributions: (i) The first term is merely the 
Coulomb interaction between the charged interface and the test ion 1. (ii) The second 
term involves the interaction between the test ion 1 and the remaining solute ions, with 
9 apif i)'f^2) being the pair distribution function and P/3{z2) the local ion density. If the 
former is approximated by ga(3{'^i,r2) ~ 1 then equation ( fT2l) becomes: 

2 



ViWaiZi) = -Vi 



(13) 



=qai'{zi) 

so that the potential of mean force reduces to the MEP times the charge, which is 
precisely the PB approximation. In other words, the PB theory neglects the (lateral) 
ion-ion correlations in that sense that 5'«/3(ri,r2) = 1^ It is for this reason that the PB 
theory is a mean-field one. Recalling that 

g^,{n.,zuz,)= ^_^,^_)^^(^^) . (14) 

where Pap{ri2, Zi, Z2) is the two-particle density function, one can equally well provide a 
geometrical interpretation: The probability of finding two ions anywhere in the solution 
is independent of their relative separation in the PB frameworkjj 

2.1.2. Debye Hiickel theory In general the strongly non-linear PB equation (jl]) can not 
be solved analytically, and its linearized version is therefore employed instead. This 
latter approach was historically first developed by Debye and Hiickel [6]. When the 
MEP is everywhere small (i.e., \eilj\ < 1), the PB equation reduces to 

Alp = (15) 

and the corresponding solution reads: 

Mz) = ^se-^' = _^2^e-"^ (16) 

e 

f Note that the existence of lateral ion-ion correlations [i.e., gaf3{ri,r2) — 1 7^ 0] have two physical 
origins: (i) Electrostatics and (ii) steric effects due to excluded volume. The latter are implicitly ignored 
in the PB framework. 

I Clearly, the bare Coulomb pair force between all constitutive ions are properly taken into account 
in the PB theory, see equation (fTS]) . It is the assumption of a structureless lateral arrangement of the 
ions that creates the crucial inconsistency in the PB framework. 
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where ips denotes the surface potential. This result can be obtained either by directly 
solving the DH equation f|T5l) or substituting the small ips value in the full PB solution 
equations ([5]) and (El). 

In order to more deeply understand the physical meaning of the linear 
approximation, we shall rewrite the DH equation (ITSll in an equivalent integral equation 
form. In this context, it is instructive to use an approximative closure for the (exact) 
Orstein-Zernike equation, as done by McQuarrie for a bulk electrolyte [9], which leads 
to the DH description: 



hoa{Zl) = Coa{Zl) + P/3 / hoi3{z2) Ci3a{ri2, Zl, Z2) d^r2 (17) 

f^l J ' ^ ' 

(DH) ' (DH) 

« -OLZx/b f» -al3lB/ri2 



,{Zl)^ + '^Pf^ j ^^^(f^ hf3a{ri2, Zi, Z2)d^r2, 



, (DH) 



where, the subscript "0" in equations (IT7I) and (ITS]) stands for the charged interface (that 
can be envisioned as a special particle species at infinite dilution). Thereby, this notation 
preserves nicely the analogy with the bulk case. The DH theory is readily obtained upon 
assuming that the direct correlation function is equal to the (sign reversed) reduced pair 
potential [i.e., c(r) = —[3V{r)], which becomes exact when poC\ and b/is oc. 
In practice it is the first line (ITTl) that is used in fluid theory to solve self-consistently 

(DH) 

the total correlation function hij = gij — 1 or the PMF via hij ^ —(3wij. 

For the sake of our discussion, however, it is the second line ( |T8|) that turns out to be 
instructive. Indeed we see now that in the DH theory, the term (7/30(^1, r2) is not trivially 
unity, since /^/^^(ri, r2) does not vanish in equation (fTSl) . in contrast to what happens in 
the PB situation. Hence ion-ion correlations are not neg lectedQ This might seem at first 
sight counter-intuitive since the DH theory is based on the linearization of PB equation 
which ignores lateral correlations. This being said, in the weak Coulomb regime where 
the DH theory is supposed to be valid, the deviations from the uncorrelated limit are 
then small. 



2.2. Strong Coulomb coupling 

2.2.1. Strong Coulomb coupling theories This last decade [TOt [TT] . a remarkable 
theoretical achievement has been accomplished in the other extreme limit of strong 
Coulomb coupling. More specifically, the counterion distribution near a charged planar 
wall has been predicted analytically and independently by Shklovskii [10] and Moreira 
and Netz [H] in the strong Coulomb coupling regime (i.e., the Gouy-Chapman problem 
at low temperature). A common and universal feature of these two works is that the 

t Note that MacQuarrie used the very same method [equation (fTSl) ] to determine analytically (via 
Fourier transformation) the DH potential in spherical geometry 9J. However, in the past, he was not 
aware of the relevance of lateral ion-ion correlations, and therefore did not point out this issue. 
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counterion distribution decays exponentially like exp{—z/b). These two approaches are 
going to be now briefly presented. 

• Using a field theoretic formulation applied to charged fluids [121 [IS] i Moreira and 
Netz [H] showed that at high Coulomb coupling (i.e.: for H = ^ ^ 1) the 
counterion distribution obeys the following exact and elegant limiting law: 

= exp(-^/6) (19) 

with (Ts = |cT|/e (having the dimension of the inverse of a surface) standing for the 
number of elementary charges per unit area. A plot of equation fll9p can be found 
in flgure [2l where a convenient visual comparison with the high temperature limit 
[equation ( ITTi) ] is offered. 

• Using a fully different and more intuitive approach, Shklovskii [T3] has applied 
Wigner crystal (WC) concepts [151 [IS] to the problem of soft charged matter at 
effective low temperature. Using some heuristic but physically sound arguments, 
essentially based on the simple fact that a "desorbed" counterion from the 
(triangular) WC counterion layer is correlated to the hole left behind over the Gouy- 
Chapman length b, Shklovskii [ll] obtains (up to the here important prefactor) 
the same result [equation (fTOjl ] as Netz. Interestingly, if one combines (i) the WC 
approach that provides the correct exponential decay exp(— z/6) and (ii) the contact 
theorem which imposes the prefactor, 2t:('.b(^1 [12] El then one recovers the exact 
answer (|T9l) . 

2.2.2. Overcharging and Thomson problem As long as the Coulomb coupling between 
ions is "fairly" moderate (which is the case for monovalent ions in aqueous solution), 
the PB theory [H [5l [181 [IS] ^"^^ even the DH one [20] describe astonishingly well the 
ion distribution (and hence the thermodynamical system properties) when compared to 
computer simulations [HI [221 [231 [211 [251 [21], theories going beyond the mean-fleld PB 
level [2Z], and even experiments [2H]- Nonetheless, as soon as ion- ion correlations get 
relevant, mean fleld theories such as the PB one [29] or its linearized version (as related 
above in 12. II) can not explain the experimentally observed relevant effect of overcharging 
[301 [31]. 

Naively, one would think that the stable conflguration corresponds to an exact 
neutralization of the macroion by the counterions. This intuition is only correct for the 
case where the counterions are uniformly smeared out over the surface of the colloid. 
Indeed basic electrostatics show that, for a central charge Z^e < (representing the 
macroion) and the shell of the counterions of radius a and (total) charge zi^'^^^^'^e > 0, 
the electrostatic potential energy is given by [32] 

E = ^I!l^l—^ + ^^——^^ (CGS) (20) 

I Note that 2t: IbctI = ' ^^^^ theory predicts the exact contact value as well [compare 

with equation pT|) ]. see figure [3 This is not true, however, for the DH version. 
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where the first term describes the interaction between the central ion and the charged 
shell and the second one is the electrostatic energy stored in the shell (i.e., the work 
done upon bringing the counterions from infinity to their current location r = a of the 
shell). Thereby, the criterion of stability 

and , = eVa > ^ Z^^'^^") = -Z^ (21) 



shows that the stable configuration corresponds to an exact neutralization. In reality, 
the counterions are discrete and not smeared out, and when electrostatically bound to 
the macroion's surface, they will maximize their separation such as to minimize the 
counterion-counterion repulsion. This problem turns out to be exactly the one that 
was addressed one century ago by Thomson p3J (also called the Thomson sphere or 
Thomson problem) who studied the ground state energy and structure of (classical) 
electrons confined on a sphere (model of a classical atom). The Thomson problem has 
only exact solutions for small and some magic numbers (e.g., N = 12 corresponding to 
the fuUerene structure) [3^. Nonetheless, based on Wigner crystal ideas [151 [TOl [35] , 
analytical model was developed which quantitatively accounts for the energy gain upon 
adsorbing overcharging counterions [36l [32]. More precisely, the following relation 
for the energy variation (relative to the globally neutral state characterized by 

n = overcharging counterion and = Zm/Zc counterions, see figure [3] for a typical 
counterion arrangement) as a function of the number n of (excess) overcharging Z^-valent 
counterions [.32J was derived: 

AE„ = -^£L \{N + nf/^ - iV3/21 + (CGS) (22) 

V47ra2 2a 

where a (~ 2) is a numerical geometrical prefactor that was determined by simulations 
(deduced from the value of AEi). [j| The first and attractive term in equation (|22ll 
stems basically from the interaction between a counterion and its oppositely charged 
Wigner-Seitz cell. Energy profiles of equation f[22|) are sketched in figure [3l where one 
can see that these analytical predictions are pretty robust. This simple approach to the 
understanding of the overcharging via the Thomson problem, Wigner crystal concept 
and computer simulations has triggered a new interest in the community [371 [3H1 [39l [ID] 
for the Thomson problem applied to soft matter. 

We now consider the problem of a pair of macroions. In [36], it was shown that 
two equally charged spheres are likely to be overcharged and undercharged in the 
strong Coulomb coupling regime leading to a metastable ionized state that yields a 
strong long-ranged attraction due to a monopolar contribution. All the mechanisms, 

f To achieve overcharging in nature one should normally add salt to the system to ensure global 
electroneurality. For the sake of simplicity, however, we will consider non-neutral systems because they 
can, on a very simple basis explain, why colloids prefer to be overcharged. 

I Note that in the case of vanishing curvature (i.e., a/dc — > oo where dc is the mean distance 
between counterions) our expression becomes exact since the planar WC limit is recovered for which 
a = 1.960516... M- 
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Figure 3. (a) Snapshot of the ground state structure with Zm = 180 corresponding 
to = 90 counterions. Note the local triangular arrangement on the "Thomson 
sphere", (b) Electrostatic ground state energy (in units of fc^To with Tq being the 
room temperature) as a function of the number of overcharging counterions n for three 
different bare charges Z,„. The neutral case was chosen as the potential energy origin, 
and the curves were produced using the theory of equation (|22p . compare text. Data 
taken from [5^ . 



SO far reported in the literature, can only explain short-ranged like-charge attraction 

To further rationalize this phenomenon and the stability of ionized states [521 [32] , 
two charged spheres of same radius a, carrying the same electric sign of charge but 
characterized by a charge ratio pz such that < pz = Zb/Za < 1, were considered. 
Starting from a macroion pair where each macroion is neutralized by its counterions, the 
process where a counterion is transfered from macroion B (low bare charge) to macroion 
A (high bare charge) was investigated [521 [32]. Having demonstrated that the ability of 
a macroion to get overcharged increases with growing (bare) surface charge density a 
(or the bare charge at fixed radius), it is clear that this counterion-transfer process will 
be energetically favorable below a certain value of pz- This theoretical prediction shows 
that the criterion for stable ionized states (latter also called by other authors [391 EH] 
"auto-ionization" ) is governed by the value of 

V^-V^>1 (23) 

(with Na/b = Za/b/Zc being the number of counterions of macroion A/B) which reflects 
the correlation- hole energy difference between the two macroions (at identical radii). In 
particular, it was demonstrated that the higher the charge-asymmetry (i.e., pz) the more 
stable the ionized state and concomitantly the higher the degree of ionization [521 [32] . 
The main findings related to this work [361 [521 [32] . can be summarized as follows: 

• The ground state of a charged sphere is always overcharged due to counterion 
correlations. 

• At finite temperature and in the strong Coulomb regime (accessible with multivalent 
aqueous ions), colloids having different bare surface charge density auto-ionize due 
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Figure 4. Computer simulation snapshots of counterion ground state configurations. 
The discrete colloidal surface charges are in white. The counterions are in blue, (a) 
"Low" and (b) "high" surface charge densities are shown. Data taken from [55] . 



to counterion correlations. 
2.3. Discretely charged surfaces 

The structural (i.e., bare) charge of spherical macroions is usually modeled by a central 
charge, which, by virtue of the Gauss' law, is equivalent to a uniformly charged macroion 
surface as far as the electrostatic field (or potential) outside the sphere is concerned. 
However, in nature the charges on the colloidal surface are discrete (exactly as the 
counterions are) and localized, see figure HI Thus, a natural question that rises is: Why 
and how does the counterion distribution depend on the way the structural charge of 
the macroion is represented (i.e., uniformly charged or discrete charges on its surface)? 
It is precisely this problem that was addressed in [SU [55] . 

Why is the counterion distribution sensitive to the choice of the representation of 
the macroion charge (discrete vs. uniform)? This question can be best answered by 
looking at and comparing the (intrinsic) electrostatic potentials generated by discretely 
and uniformly charged macroions (without counterions) [M]- It was demonstrated in 
Ref. [51] that the electrostatic potential at a reduced distance r/a from the sphere 
(where a stands for the distance of closest approach between an external unit test-charge 
and the macroion surface) may be significantly different according to the nature of the 
macroion charge. In particular we show that the higher the bare surface charge (i.e., 
the closer we get to a uniform charge distribution) the shorter the correlation length 
(typically ~ a/T/o^) between the discrete surface charges, as intuitively expected. 
More specifically, the contact potential is sensitive to the localization of the discrete 
charges, leading to a pronounced depth in their vicinity. All those features, solely based 
on the spatial behavior of the electrostatic potential stemming from the bare macroion, 
indicate that the counterion distribution should be much more complicated for a discrete 
macroion surface charge distribution than for the uniform case. 

We now come to the other important question: How is the counterion distribution 
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modified when introducing the more reahstic discrete macroion's surface charge 
distribution? This point is thoroughly addressed in |55], where two regimes are 
considered: Ground state (T = 0) and finite temperatures. The corresponding relevant 
findings [SH [55] can be summed up as follows: 

• At zero temperature, the counterion (surface) structure possesses greater order the 
the higher the reduced surface charge density as and/or counterion valence Zc are. 

• When overcharging comes into play several scenarios occur: (i) At large as, the 
overcharging is nearly the same as that obtained at a uniformly charged macroion's 
surface, (ii) At low ag and for monovalent counterions, overcharging is always 
weaker for discrete macroion charge distribution, due to the ion-pairing frustration 
for the excess counterions. (iii) At low as and for highly multivalent counterions, 
overcharging can even be stronger in the discrete case due to ion-pairing. 

• At finite temperature (in aqueous solutions), the volume counterion distribution is 
only affected for low as and multivalent counterions. 

The effect of surface charge discretization was later examined for different 
geometries by several groups [56l ETJ EH [59], [60] [61] [62] . 

3. The crucial role of excluded volume 

3.1. Monovalent ions near a charged sphere 

So far, we have had a pretty good understanding of the physics involved in the counterion 
distribution for salt-free systems where excluded volume effects are irrelevant. The 
situation becomes much more complicated at finite salt-concentrations in aqueous 
solutions (i.e., water at room temperature in the presence of added salt), where the 
Coulomb coupling is (rather) weak especially for monovalent ions. Thereby, a direct 
application of Wigner crystal ideas is not straightforward enough to account for the 
unexpected overcharging at weak Coulomb coupling that was reported theoretically 
[631 [Ml [651 [66]). but unexplained, for monovalent salt-ions of large size. 

Molecular dynamics computer simulations as well as integral-equation theory 
[67] were employed to identify the mechanisms that govern counterions ordering and 
overcharging in this weak Coulomb coupling regime. Those mechanisms are as follows: 

• Increasing the electrolyte particle size (at given salt concentration) decreases the 
available volume of the fluid (or equivalently its entropy) which favors ion-ion 
correlations. 

• The interface provided by the macroion causes an increase of the ion density close 
to it, and concomitantly enhances the lateral ordering (similar to the prefreezing 
phenomenon in neutral inhomogeneous fluids). 

• Surface lateral ordering and (weak) Coulomb coupling lead to overcharging. 
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Figure 5. Schematic view of the electrostatic model for macroions near an oppositely 
charged interface. The macroions are characterized by a distance of closest approach 
z = a to the charged surface, leading to two screening strengths kq and k for < z < a 
and z > a, respectively. 



3.2. Macroion adsorption at planar substrates 

Excluded volume effects coupled to electrostatic interactions can also lead to counter- 
intuitive phenomena in the process of macroion adsorption. A description of the model 
setup is sketched in figure El For instance, Jimenez- Angeles and Lozada-Cassou showed, 
using integral equation theory, that for moderately (attractive) charged substrates, a film 
of coions first builds up. The electrostatic consequence is that at the direct vicinity of 
the surface of the substrate its charge gets amplified (i.e., surface charge amplification). 
The driving force of this effect is due to the macroion-ion attractive correlations. This 
effect was overlooked in the past, because the authors neglected either the finite size of 
the macroion [68] or the spatial distribution of the little salt-ions [69] . 

Recently, this problem was revisited by using a very simple analytical model 
based on the Debye-Hiickel approximation hut taking into account the finite size of 
the macroion via its distance of closest approach a (i.e., its radius) to the wall (see 
figure [5] as well) [70|. Two regimes were specifically examined: The strong and weak 
screening regimes which are now briefly described. 

• In the strong screening regime {kqO ^1)0 the wall-macroion attractive interaction 

f The negative counterions of the positively charged macroions correspond to the coions of the planar 
substrate. Thereby, electrostatic correlations tend to localize the counterions of the macroions over its 
whole surface in a uniform manner. Hence, as long as the strength of the surface charge density of the 
oppositely charged substrate is low enough, a finite number of counterions of the macroions should stay 
in the vicinity of the interface (see figure [5]), leading to a surface charge amplification, 
f stands for the screening strength stemming uniquely from the little ions, see also figure [H 
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is exclusively governed by the screening contrast kq/k. [|| More precisely, it was 
shown that the contact potential of interaction Um is merely given by [70] : 



PUm^l- -■ (24) 

In the weak screening regime {koO <^ 1) and for sufficiently small surface charge 
density ^ 1), the reduced electric field [§| at contact follows this simple law 

This equation fl25l) tells us that surface charge amplification is increasing with 
growing colloidal particle size a and increasing Gouy-Chapman length b (i.e., 
decreasing as). 



4. Image charges in spherical geometry 

In a typical experimental setup, the dielectric constant of a macroion is rather low 
{Em ~ 2 — 5) which is much smaller than that of its embedding solvent (e. g., for water 
^soiv ~ 80) leading to a high dielectric contrast, = at the interface. It turns 

out that for a perfect planar substrate (which can be envisioned as a colloid of vanishing 
curvature), there is an elegant analytical solution for the electric field. More precisely, 
the electric field generated by the induced surface charge at the interface positioned at 
z = (due to the presence of a point-like ion of charge q located at z = i) can be exactly 
obtained by a "fictive" point-like charge = A^g located at the mirror position z = —i 
\71\ . This feature corresponds to the so-called method of image charges. The inclusion 
of such image forces for the case of an electrolyte close to a planar dielectric interface was 
studied in the past by computer simulations [72], [731 [711 [56] , integral equation formahsms 
[75| [76] . mean- field [771 EH UM HSl EH! and strong-coupling [56] theories. As far as the 
cylindrical case [HI [82l [83] is concerned, there is no simple "image charge" picture. 

The problem of the dielectric discontinuity in spherical geometry is, already at the 
level of a single ion interacting with a dielectric (neutral) sphere, considerably more 
complicated than its planar counterpart. Indeed, if we want to reformulate the problem 
in terms of image charges, one would need an infinite number of image charges, thus 
making its usage much less attractive than in the planar case. Due to this difficulty, 
the problem of image charges in spherical geometry is sparsely studied in soft matter. 
Nevertheless, twenty years ago, Linse studied the counterion distribution with image 
forces around spherical charged micelles by means of Monte Carlo simulations [HI]. In 
his work [84j, Linse used a two-image charge approximation instead of the full continuous 
image charge distribution. The conclusions of his study remain qualitatively correct. 

I K stands for the total screening strength stemming from aU the ions present in the solution (also 
including the macroions), see also figure [5l 

§ The reduced electric field is defined as E*{z) = — (6/2)e-^^^^ such that at the interface z = we 
have E*{0) = -1. 
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Figure 6. (a) Model for a dielectric sphere (colloid) of dielectric constant Sm 
embedded in an infinite medium (the solvent) characterized by a different dielectric 
constant Ssoiv A test positive charge (q) is located near the boundary outside the 
spherical particle at a radial distance R. The resulting induced surface polarization 
charges are also illustrated for the case where Ssoiv > £m ■ Note that the global induced 
net charge vanishes. This is a two-dimensional representation of the three-dimensional 
system, (b) Polar profiles, as obtained from equation (pS]) . of the surface density of 
polarization charge Upfi'^^O) in units of cr^°] — ^^^^^ ^2 for different radial distances 
R of the test charge q with Ssoiv — 80, = 2 and a — 7.5c?. 



The dielectric response of a dipolar fluid confined to a spherical cavity was recently 
addressed by Blaak and Hansen using MD simulations [85] . 

In the field of image forces in spherical geometry, exact results for the electrostatics 
of an ion interacting with a dielectric sphere (see figure [6] for the model geometry) were 
reported [86]. Furthermore, Monte Carlo simulations were performed to elucidate the 
behavior of an electrolyte near a spherical macroion at finite dielectric contrast, where 
image forces are properly taken into account [86]. The main results are as follows |86j : 

• Single ion: A compact and exact analytical expression has been derived for the 
polar profile of the induced surface charge, and it reads: 

AtTEsoIvR^ ^ esolv[l + 1) + 

where g is a test ion at a radial distance R (see figure [6]) and Pi designates the 
Legendre polynomials of order I. The strength as well as the range of image forces 
in spherical geometry are considerably smaller than at vanishing curvature, due to 
auto -screening. 

• Electrolyte: For aqueous monovalent ions the (effective) image force is basically 
equal to the self-image one (i.e., the interaction between an ion and its own image). 
However, when dealing with multivalent counterions, the lateral image-counterion 
correlations can significantly affect the (local) counterion density and, as a major 
effect, they screen the self-image repulsion. Upon adding salt, it was shown that 
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Figure 7. A computer simulation snapshot of PE-colloid complexation (tennis ball- 
like conformation) [5T]. 

the strength of the image forces induced by the coions is marginal. Besides, 
overcharging is robust against image forces. 

Very recently, Rescic and Linse extended [87] the one-colloid problem to the two- 
colloid interaction problem with dielectric discontinuity. Using a cylindrical cell model 
and MC simulations, they found (i) weaker counterion accumulation at the macroion's 
surfaces, (ii) stronger effective repulsion at moderate Coulomb couphng, and (iii) a less 
attractive effective force at strong Coulomb coupling. These findings are fully consistent 
with the one-colloid features just discussed above. 

5. Polyelectrolyte adsorption and multilayers 

Polyelectrolytes (PEs) are polymers containing a variable (usually large) amount of 
ionizable monomer along the chemical backbone. Once dissolved in a suitable polar 
solvent such as water, the ion pairs dissociate by creating a charged chain with floating 
counterions. PEs represent a broad and interesting class of materials that have attracted 
an increasing attention in the scientific community. PEs have applications in modern 
technology as well as biology, since virtually all proteins, as well as DNA, are charged. 
The adsorption of PEs onto surfaces is an important process, since they modify the 
physico-chemical properties of the surface. From a theoretical point of view, charged 
polymers (in bulk or adsorbed) are much less understood [8811H9] than neutral ones [90] . 
One of the main difficulties is the addition of new length scales set by the tremendous 
long-ranged Coulomb interaction. Hence, the study of adsorption of PEs is motivated 
by fundamental aspects as well as practical ones. 
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5.1. Poly electrolyte- colloid complexation 

5.1.1. Oppositely charged spherical substrates The works related to the interaction 
between PEs and oppositely charged spheres are here briefly reviewedl^l The 
complexation of flexible PEs with oppositely charged macroions is a relevant process 
in biology [9l]. For instance a nucleosome can be seen as an electrostatic binding 
between DNA and histone proteins, where the latter can be envisioned as charged 
spheresfl Many theoreticians [951 EHl [971 EHl EHl M [Ml [M [921 M [M [TnH] 
have investigated these types of objects to understand the electrostatics governing 
those structures. Two very relevant results are: (i) the possible overcharging of the 
sphere by the long PE and (ii) a strong wrapping of the PE about the sphere (see 
figure [7] for an example). A considerable effort was also provided by the simulators 
[Tn6t flOTl [TM fTn9l fTTOl [TTTl fTT2l [ml fTM fTT5l fTTHl M, W\ these last twelve years. 
Some relevant findings in this field can be summarized as follows: 

• The effect of chain stiffness, which was first systematically studied by Wallin and 
Linse |106j by MC simulations, is an important key controlling the PE adsorption. 
They showed that the lower the chain stiffness, the higher the PE adsorption 
and, concomitantly, the overcharging of the charged sphere by the PE. Stoll and 
Chodanowski |114j . using MC simulation as well but with Yukawa potentials, 
showed that upon increasing the chain stiffness, solenoid conformations are obtained 
as predicted analytically by Nguyen and Shklovskii [lOlJ. 

• The effect of chain length was also addressed (by means of computer simulations) 
in the past |108[ I112[ 11131 1115] . For large chain/sphere size ratio, Chodanowski 
and Stoll |112j found, for fully flexible chains, that only a marginal portion of the 
PE gets adsorbed to the sphere, and the rest of the chain consists of extended 
tails. At "moderate" chain/sphere size ratio [112j . they found a strong PE collapse 
into a tennis-ball like structure (as illustrated in figure [Tj). Considering both the 
effects of chain length and the chain stiffness, Akinchina and Linse |115j reported a 
rich phase behavior: Tennis-ball like, solenoid, multiloop (also called rosette |94j). 
single loop, as well as "U"-shaped conformations. Note, that there is remarkable 
agreement with the rosette structure found theoretically by Schiessel et al. |118] . 

• The effect of the discrete nature of the protein charge distribution was addressed by 
Carlsson et al. [119] . In their MC simulations 1119], they found that complexation 
can be stronger with a discrete protein charge distribution (in agreement with the 
ideas discussed in 12.31) . 

• Multisphere complexation involving many charged spheres bridged via oppositely 
charged PEs were investigated by Jonsson and Linse |110[ 1111] by means of MC 
simulations. The effect of linear chain charge density, chain length, and macroion 



f The reader can also look at recent reviews [92l [93] on this field for more details. 
I We are aware that this assumption is at best a caricature of a real system (provided that non-specific 
interactions are dominant). Nonetheless, from an electrostatic viewpoint, we think that the qualitative 
features should be captured. 
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charge valency was addressed in Ref. |llUj . Interestingly, at prescribed PE linear 
charge density, the authors found that complexation gets stronger upon increasing 
the chain length [llOj . The effect of chain flexibility was studied in jlllj , and it was 
found that the macroion arrangement gradually becomes more linear and ordered 
along the (long) chain when its stiffness is increased. 

5.1.2. Like-charge complexation Whereas many studies have been devoted for the case 
of chain-sphere complexation where the two charged bodies are oppositely charged, 
as we just saw, much less is known concerning the problem of like-charge sphere- PE 
complexation. 

In |120l 1121] . the complexation between a sphere and a long flexible PE {both 
negatively charged) was discussed. Whereas like-charge attraction in the strong Coulomb 
coupling limit is expected (and therefore complexation too), new and rather unexpected 
chain conformations are reported. Different coupling regimes as well as the influence 
of the linear charge density, /, of the PE chain were considered in |121j . The relevant 
conclusions are as follows: 

• At strong coupling the PE chain is always adsorbed in a flat structure, whose 
conformation strongly depends on /. At high /, the conformation consists of 
densely packed monomers following a Hamiltonian-walk. Upon reducing / the 
chain tends to spread more and more over the particle surface. These findings 
could have some relevance for organic solutions. 

• Under aqueous conditions, complexation can be obtained with multivalent 
counterions and for high enough values of /. In contrast to the strong coupling 
case, the formation of loops is reported. 

5.2. Poly electrolyte adsorption at planar surfaces 

The reader who wants to know a detailed account of the field of PE adsorption at surfaces 
is invited to consult the recent reviews of Netz and Andelamn |122j and of Dobrynin and 
Rubinstein [123] . In this part, one would like to propose some basic ideas and features 
supported by MC simulations about the adsorption of highly charged polyelectrolytes 
onto oppositely charged planar surfaces in a salt-free environment [1241 11251 I126j . 
Flexible |124[ 1125] as well as rod-like |126] PEs are now discussed. 

5.2.1. Role of entropy There is a simple and clear entropic mechanism that influences 
multi-polymer- chain adsorption that is going to be pointed out first. It can be best 
understood by recalling the counterion release effect: The adsorption process of one 
polyion of valence Z typically leads to the release into solution of Z (initially adsorbed) 
surface monovalent counterions, which is "electrostatically invariant" but entropically 
(highly) favorable. This very same effect is also the reason why longer chains can 
better adsorb at a prescribed monomer density. Indeed, at prescribed monomer density. 
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increasing the chain length ivj^ involves decreasing the the number of chains. Thereby, 
the resulting (bulk) entropy stemming from the PE chains becomes reduced accordingly. 
This entropic mechanism linked to the chain length at prescribed monomer density is 
henceforth referred to as: polymerization induced adsorption. 

5.2.2. Flexible chains [124] 

When no image forces are present (i.e., = 0), it was found that the monomer 
density profile, n{z)., decays monotonically for very short chains even near contact, see 
figure [H](a). Longer chains experience a short-ranged repulsion in the vicinity of the 
charged wall {z < d) due to chain- entropy effects, 

When image forces come into play, (partial) monomer desorption sets in, whose 
strength increases with growing chain length Nm- This feature is due to the repulsive 
image-chain interaction that scales like A^^, whereas the attractive Wigner crystal 
correlation^ scales only like Nm"^ . 

The fraction of charge (J*{z) of the fluid as a function of monomer- wall separation, 
2;, is another interesting quantity to characterize the adsorption behavior. At = 0, 
overcharging [as signaled by cr*(z) > 1] occurs as soon as chains are longer than dimers, 
see figure 4(a) in [124j . In the presence of image forces, the strength of the overcharging 
is nearly identical to that obtained without image forces at A^ = (compare with figure 
4(a) in Ref. |125j ). Thereby, the main effect of image charges is (i) to decrease the 
fraction of charge (J*{z) near contact {z < 1.2a) upon growing Nm and (ii) to (slightly) 
shift the position of the maximum of <J*{z) to larger z. 

5.2.3. Rigid chains [T26] 

Dimers exhibit a monotonic behavior for n{z) that is similar to point-like ions. 
For longer chains there exits a small monomer depletion near the charged wall for 
an intermediate regime of Nm-, see figure [8]^a). At high enough Nm-, n{z) reveals 
again a monotonic behavior, see figure 1(a) in Ref. |126] . This interesting effect is 
the result of two antagonistic entropy-driving forces, namely, (i) chain-entropy and 
(ii) polymerization induced adsorption. Electrostatic chain-chain correlations, whose 
strength grows in a non-trivial way with A^mili] favor also chain adsorption. Figure [8]^a) 
clearly shows that the adsorption of rigid PEs is much stronger than that of flexible 

f Rigorously, Nm represents the number of monomers per chain corresponding experimentaUy to the 
polymerization degree. 

I The chain-entropy effect here is merely due to the much lower number of available conformations in 
the adsorbed state. It has to be distinguished from that previously discussed in 15.2. II 

§ When charged polymers are adsorbed on the surface, they also tend to build a Wigner crystal due 
to the strong mutual Coulomb inter-chain repulsion. The higher the chain length Nm (i.e., the length 
of the chain) the stronger the effect. At prescribed reduced surface charge density CTs and monomer 
concentration, this leads to a 2D plasma term (i.e., interchain repulsion reduced by thermal energy) 

3 /2 . 

that roughly varies like Nm , as is the case for point-like multivalent ions. 

II Due to the strong extension of the chain, it is no longer suitable to use the point-like and/or spherical 

3/2 

polyion picture leading to the WC term in Nm ■ 
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Figure 8. (a) Profiles of the monomer density ti{z) for various chain length N„i: 
flexible vs rigid chains [126j . Snapshots at N„i = 8 for (b) flexible and (c) rigid chains. 



ones. This feature is also detectable in the snapshots, see figure [8](b) andlUj^c). 

Upon polarizing the interface, it is found that the degree of adsorption is 
considerably reduced. Nonetheless, a comparison with the flexible case |125] shows 
that the values at contact at flnite are quite similar. 

5.2.4- Summary To sum up, MC simulations [1241 11251 1126] show us that: 

• Without a dielectric discontinuity (A^ = 0), flexible PE chains experience short- 
ranged repulsion near the charged substrate due to chain-entropy effects. In 
contrast, rigid PE chains are more strongly adsorbed (due to a weaker loss of 
chain-entropy) and, when long enough, experience a purely effective attraction. 

• Image forces lower the degree of adsorption for flexible and rigid PE chains. 
However, the overcharging of the substrate by the PEs is robust (irrespective of 
the chain flexibility) against image forces. 

5.3. Poly electrolyte multilayering 

PE multilayer thin fllms are often obtained using a so-called layer-by-layer deposition 
technique |127l I128j : A (say negatively) charged substrate is alternatively exposed to 
a polycation (PC) solution and a polyanion (PA) one. This method and the resulting 
materials have a fantastic potential of application in technology, e. g., biosensing |129j . 
catalysis [130] . nonlinear optical devices (131j . nanoparticle coating [132] . etc. 

From the theoretical side the literature is rather poor. However, a few analytical 
works about PE multilayers on charged planar surfaces based on different levels of 
approximation are available |133l 11341 1135j . Solis and Olvera de la Cruz considered the 
conditions under which the spontaneous formation of polyelectrolyte layered structures 
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can be induced by a charged wall |133j . Based on Debye-Hiickel approximations for the 
electrostatic interactions, but including some lateral correlations by the consideration of 
given adsorbed PE structures, Netz and Joanny [l34] found a remarkable stability of the 
(semi-fiexible) PE multilayers supported by scaling laws. For weakly charged flexible 
polyelectrolytes at high ionic strength, qualitative agreements between theory [135j . 
also based on scaling laws, and experimental observations |136j (such as the predicted 
thickness and net charge of the PE multilayer) were achieved. More recently, Shafir and 
Andelman, using mean-field theory, pointed out the relevant role of a specific strong 
short-range interaction between PAs and PCs. 

A tremendous difficulty in PE multilayering is the strong electrostatic correlations 
between PCs and PAs, which are hard to be satisfactorily taken into account in 
(modified) mean- field theories. In this respect, computer simulations are of great 
help. The first simulation model for PE multilayering was developed in [91]. Later 
Panchagnula et al. performed similar computer simulations [137] . where the dynamical 
aspect was more emphasized. Several types of substrate geometry were considered, from 
spherical particles [911 1137[ I138j to planar substrates |139i 1140] via cylindrical ones [141] . 
Relevant simulation findings for spherical [9lj and planar substrates |139j are going to 
be described. 

5.3.1. Poly electrolyte multilayering at spherical substrates From the study in [91] 
concerning substrates with finite radii (i. e., charged spheres), one has learned that 
non- electro static forces are required to obtain (true) PE multilayers. More precisely, by 
introducing a (additional) short-range van der Waals-like attraction (whose strength is 
characterized by its value at contact, Xvdw, in units of fc^T) between the substrate's 
surface and the (monomers of the) oppositely charged chains. The PE structure results 
then from a complicated interplay between: (i) PC-PA strong attraction (favoring a 
collapse into a compact globular state) and (ii) PE-substrate correlations (favoring flat 
adsorption and wrapping^^ around the sphere). Briefiy, the main findings in [91] are as 
follows: 

• Flat 62/a?/er-structures, involving two long oppositely charged chains, set in only 
for large enough Xvdw At low Xvdw, the strong driving PA-PC force leads to PE 
globular structures, see figure [91 

• Stable and fiat multilayers are only obtainable at large enough Xvdw In a purely 
electrostatic regime {Xvdw = 0) PE globules are formed preventing a uniform 
coverage of the surface, see figure [lOl 

• Short chains are not suitable candidates for PE multilayers, due to (i) the weaker 
effect of polymerization adsorption and (ii) reduced chain-chain correlations. 

f Note that the wrappmg from the chain(s) around the colloid is peculiar to spherical substrates. 
Besides it should be reminded that wrapping is also governed by the repulsive interaction between the 
turns of a chain [351. 
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Figure 9. Typical configurations for one PC (in white) and one PA (in red) adsorbed 
onto the negatively charged colloid at different Xudto-couplings. (a) Xvdw = (b) 
Xvdw = 1 (c) Xvdw = 2 (d) Xvdw = 3 (e) Xvdw = 5. Note the remarkable structural 
change occurring at Xvdw = 2. The small univalent counterions (anions and cations) 
are omitted for clarity. 




Figure 10. Typical equilibrium configurations for 12 PEs (6 PCs in white and 6 
PAs in red) adsorbed onto the negatively charged colloid at different Xi)du)"COuplings. 
The little counterions (anions and cations) are omitted for clarity, (a) Xvdw = (b) 

Xvdw — 3. 



5.3.2. Planar substrates PE multilayering onto planar substrates were investigated in 
|139[ 1142] 1140] . The zero- curvature case differs qualitatively from the spherical one. First 
the intrinsic electric field is higher in the former case 0. Secondly the chain- wrapping is 
no-longer present at zero curvature. Consequently at given surface charge density, we 
expect a stronger PE-layering. The important results can be formulated as follows: 

f At zero curvature we have -0 ~ r in contrast to finite curvature where ■0 ~ ^/r 
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Figure 11. Profiles of monomer density n±{z) for oppositely charged PEs adsorbed 
onto a negatively charged planar substrate. Xvdw-covcp\ings. (a) Xvdw = (b) Xvdw = 5 
|139j . The snapshots shown correspond to chain length Nm = 20. 



• As with for spherical substrates, the relevance of short-ranged non-electrostatic 
forces is also demonstrated here, see figure [TTJ Flat multilayers can not be achieved 
with solely electrostatic forces. 

• The formation of islands (i.e., clusters of PC-PA chains) onto the substrate are 
reported [139j and qualitatively confirm the experimental observations of the early 
stages of PE deposition (one or two bilayers) |143l I144j . 

6. Confined crystalline colloids 

It is well known from solid state studies that strongly confined (i.e., quasi two- 
dimensional) systems exhibit properties and a phase behavior that may drastically 
differ from those in the bulk [145] . This feature is also vivid in colloidal systems, 
and those materials represent ideal model systems to analyze (experimentally as well as 
theoretically) and understand confinement effects on a mesoscopic scale corresponding 
to the interparticle distance. Using external fields, a colloidal system can be prepared 
in a controlled way into prescribed equilibrium and non-equilibrium states |146j . For 
instance, in equilibrium, sohdification near interfaces (provided by a substrate or a large 
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"impurity" ) can occur under thermodynamic conditions where the bulk is still fluid (so- 
called prefreezing) . In non-equilibrium, a wall may act as a center of heterogeneous 
nucleation (favored by the excess surface-energy already offered by the wall/nucleus 
interface) and initiate crystal growth. Most of our experimental knowledge of freezing 
in a conflning slit-like geometry is based on real-space measurements of mesoscopic 
model systems such as charged colloidal suspensions between glass plates |147[ 1148] . 

In this section, different relevant achievements in the fleld of confined charged 
colloidal crystals are discussed. 

6.1. Two-dimensional dipolar mixtures 

Two-dimensional colloidal systems can be achieved for instance via sedimentation and 
trapping at the air/water interface |149[ 1150] . At strong external field, all the dipolar 
moments are aligned in the direction of the applied external field, leading to a purely 
repulsive pair interaction that scales like: 



where mi and m2 stand for the induced dipolar moments of the particles 1 and 2, 
respectively. 

Whereas the one-component situation trivially yields a triangular lattice, the 
binary mixture situation provides a very rich phase behavior [152] . This feature can 
be conveniently exploited for potential technological applications: optical band-gap 
materials (so-called photonic crystals) [153] . molecular sieves |154j . nano-filters with 
prescribed porosity |155j . etc. There have been recent advances in this field that are 
going to be concisely explained here. 

Two dimensional binary mixtures made up of two types of dipolar particles [(i) big 
particles with a large dipolar moment (species A) and (ii) small particles with a small 
dipolar moment (species B)] were investigated experimentally [150] . The corresponding 
setup and a representative snapshot of the microstructure are shown in figure [121 A 
remarkable feature is the stability of the square phase at strong dipolar asymmetry 
{mB/niA ~ 10%). 

On the theoretical side, the phase behavior of such a binary dipolar mixture 
at zero temperature was studied using lattice sums [ 152] . The relevant reduced 
parameters are (i) the reduced dipolar moment m = mB/mA and (ii) the composition 
^ = ns/iuA + ub). The resulting phase diagram is shown in figure [131 The main 
results are as follows: 

• The phase diagram qualitatively differs from that of hard disks [156j . For low 
dipolar asymmetry m > 0.5 a stable mixture T{AB2) sets in (see figure [T3|) in 

f Note that in the experimental situations, one has often to deal with magnetic colloidal particles (so 
called ferrofluids). However electric dipoles are also realizable '1511. This being said, regardless of the 
nature of dipolar moment (i.e., magnetic or electric), the same pair interaction (|27p prevails. Hence 
results on (super)magnetic particles enter also adequately in the scope of this review. 
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Figure 12. (a) Super-paramagnetic colloidal particles confined at a water/air 
interface due to gravity. An external magnetic field H perpendicular to the interface 
induces a magnetic moment m in each bead leading to a repulsive dipolar interaction, 
see equation ([27)) . (b) Micrograph showing three touching square-latticed grains at low 
reduced temperature with a global composition X = ns/inB + nA) ~ 45% [with ua 
and Ub standing for the area density of the big and small particles, respectively] and 
a reduced moment mB/mA ~ 10%. Data taken from |150j . 
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Figure 13. The phase diagram in the (rn,X)-plane at T = 0. Three important 
phases are shown: S{AB), T(A2-B) and T{AB2). The reader can find more details 
about the other structures in |152j . Data taken from [152j . 



contrast to the case of hard-disk mixtures where no mixture is predicted at low size 
asymmetry |156] . The stabihty of this phase T{AB2) was also reported in molecular 
dynamics simulations [I57j . At even smaller dipolar asymmetry m > 0.88, an 
additional (globally triangular) phase mixture T{A2B) is stable, see figure [131 

• The stabihty of the square phase S{AB) (see figure [13]) is in good qualitative 
agreement with the experimental findings in |150j . where the dominance of the 
square phase is also reported (see figure [12]) as previously mentioned. 
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6.2. Crystalline colloidal bilayers 

Crystalline bilayers made up of charged particles have been intensively studied these 
last years in the soft matter colloid community [1581 1159] as well as in the solid state 
physics (classical |160l 11611 11621 11631 1164] and non-classical electrons |165l 11661 1167] ) 
and dusty plasma communities [1681 1169] . 

The effective interaction between these constitutive mesoscopic macroions is neither 
hard-sphere like nor purely Coulombic, but it is rather described by an intermediate 
screened Coulomb [also called Yukawa or DLVO (Derjaguin-Landau-Verwey-Overbeek) 
|170[ 1171] ] due to the screening mediated by the additional microions present in the 
system. The screening strength can be tuned by varying the microion concentration: 
For colloidal systems, salt ions can be conveniently added to the aqueous suspension; 
The dusty plasma, on the other hand, consists of electrons and impurity ions. 

6.2.1. Equilibrium The equilibrium phase diagram at zero temperature of crystalline 
bilayers was investigated theoretically in [158]. The constitutive (point-like) particles 
interact via a Yukawa pair potential of the form 

Vy^,{r) = Vo^-^^^, (28) 
nr 

where Vq sets the energy scaleqj The choice of this potential is motivated by the 
experimental model systems described above. The crystalline bilayer consists of two 
(identical) layers containing in total particles in the [x, y) plane. The corresponding 
(total) surface density p is then given by N/A, with A being the (macroscopic) layer 
area. The distance D, separating the two layers in the 2;-direction, is prescribed by an 
(implicit) external potential confining the system. 

The zero-temperature phase behavior is fully determined by two dimensionless 
parameters, namely the reduced layer density, rj = pD'^/2, and the reduced screening 
strength, A = kD. Using a straightforward lattice sum technique, the phase diagram 
was calculated for arbitrary A and rj, see figure finl^ The most interesting findings [158] 
are as follows: 

• Whereas the two known extreme limits of zero |163[ I162[ 1160] and infinite 
|172[ 11731 1174] screening strength A are recovered by lattice sum calculations [158] . 

f Note that in the ground state, i.e. at rigorously zero temperature, the value of Vq is irrelevant. 
Nonetheless in experimental situations, the energy amplitude Vq = ^^'^ e{i+KK)'^ scales like the square 
of the charge Z of the particles with a physical hard core radius R reduced by the dielectric constant 
e of the solvent {e — 1 for the dusty plasma). For a charged colloids, Z is typically of the order of 
100 — 100 000 elementary charges such that V{r = d) can be much larger than /csT at interparticle 
distance (d), justifying formally our zero-temperature calculations. 

I Note that the ground state at vanishing screening A — > corresponds always to bilyaers. Indeed, 
two equally charged walls do not generate anay electric field within the slit, and consequently they do 
not alter the stable Wigner crystal structure obtained at any other surface charge density (including 
neutral walls). Thereby, if one considers the special case of two walls corresponding to neutralizing 
backgrounds, the ground state is always a bilayer. At finite screening A 7^ 0, however, the situation is 
more complicated, and multilayers (i.e., beyond bilayers) are stable at high enough density rj. 
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Figure 14. The hard sphere hmit A — > oo is sketched on top. The dashed (solid) 
lines denote continuous (discontinuous) transitions. The filled region corresponds to 
the coexistence domain of phases IV and V. The vertical arrow indicates the double 
reentrant behavior of phase IVA. The insets show the lattice geometries, where the 
filled (open) circles correspond to the lower (upper) layer. 



it is demonstrated that, at intermediate A, the phase behavior is strikingly different 
from a simple interpolation between these two limits. First, there is a first-order 
coexistence between two different staggered rhombic lattices (IVA and IVB in figure 
[T^ differing in their relative shift of the two unit cells. Second, the staggered 
rhombic phase IVA exhibits a novel reentrant effect for fixed density and varied 
screening length, see figure [TH Depending on the density, the reentrant transition 
can proceed via a staggered square III or a staggered triangular solid V including 
even a double reentrant transition of the rhombic phase IVA, see figure [TH 

• A comparative study |159] of the phase behavior of highly charged colloidal spheres 
in a confined wedge geometry reveals semi-quantitative agreement between theory 
and experiment. 
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6.2.2. N on- equilibrium The non- equilibrium cas^ at finite temperature as driven by 
a linear shear flow has been addressed in [1751 1176] . The steady state developed under 
shear as well as the relaxation back to equilibrium after cessation of shear were analyzed 
with the help of non-equilibrium Brownian dynamics. The pertinent results are: 

• For increasing shear rates, the following steady states are reported: First, up to a 
threshold of the shear rate, there is a static solid which is elastically sheared. Then, 
at higher shear rates the crystalline bilayer melts, and even higher shear rates lead 
to a reentrant solid stratifled in the shear direction. 

• After instantaneous cessation of shear, a nonmonotonic behavior of the typical 
relaxation time is found. In particular, application of high shear rates accelerates 
the relaxation back to equilibrium since shear-induced ordering facilitates the 
growth of the equilibrium crystal. 

• The orientation of a crystalline bilayer can be tuned at wish upon applying a 
(strong) shear rate in the desired direction and subsequently letting the system 
relax. 

7. Conclusions 

Various electrostatic effects in soft matter have been discussed. Generally speaking, 
charged systems are fascinating because they simultaneously involve short-ranged 
excluded volume effects (as soon as the latter are properly taken into account) already 
present in neutral systems, and additionally the long-ranged Coulomb interaction. The 
latter constitutes a formidable theoretical challenge. 

In terms of similarities with classical solid state physics and (elementary) quantum 
chemistry, two striking analogies were identifled: (i) The overcharging occurring at 
a sphere reduces to the old Thomson problem; (ii) The ground state of two spherical 
macroions is ionized, with the degree of ionization (and therefore the attraction strength) 
growing with the difference in surface charge density between the two macroions. This 
behavior is highly reminiscent of the (molecular) ionic bonding jl77] where the difference 
in electronegativity between the two atoms governs its stability. 

Excluded volume effects are equally important to fully understand phenomena 
like overcharging (i.e., surface charge reversal) and surface charge ampliflcation. For 
overcharging, the counterion layer can reach a high ordering when the local packing 
fraction is raised, by simply increasing the size of the adsorbed counterions. 

Image forces turn out to be systematically short-ranged. Their effects are only 
vivid close to the substrate at distances corresponding roughly to the linear size of 

f The starting unsheared configuration corresponds to a staggered square lattice with a reduced surface 
particle density 77 — 0.24 and a reduced screening strength A — 2.5. Two walls are present to ensure the 
confinement. To this end, screened Coulomb and short-ranged (of the Lennard- Jones type) repulsive 
potentials were tested, and it was found that our results are marginally sensitive to the choice of the 
repulsive wall-particle interaction. 
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the microion^ (counterions and/or charged monomers). As far as the adsorption 
of polyelectrolytes is concerned, there are two important driving forces that act 
concomitantly: (i) The polymerization-induced adsorption that works like the principle 
of counterion release (so entropy based) and (ii) purely electrostatic lateral correlations 
(reminiscent of the classical Wigner crystal). 

Confined colloidal crystals seem to be now pretty well understood up to bilayers. 
There is presently some experimental |178l I179j and simulational |180j evidence that, 
upon increasing the projected surface particle density, the transition from two-layer to 
three-layer structures involve four (and even more) layered crystalline structures. This 
is a problem that needs an urgent and clear understanding. 

On a more "material/engineering" level, multilayered structures can apparently 
also be experimentally obtained by combining oppositely charged colloids/micelles [181j . 
instead of polyelectrolytes. To explore this new field, a considerable theoretical effort 
would be needed to identify the parameters phase space (such as salt concentration, 
charges of the colloids and the substrates, particle size etc.) allowing the onset of such 
structures without strong clustering. 
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